Lecture 22 

Four Squares Theorem 



Pell-Brahmagupta Equation (continued) - x 2 — dy 2 = 1, □ ^ d e N, if (x, y) and 
(z, w) solutions then x<z=>y<wii and only if x + \fdy < z + Vdw. 

Theorem 78. If(x\,y\) is the least positive solution of 'x 2 —dy 2 = IwhereD ^deN, 
then all positive solutions are given by (x n , y n ) where x n + \fdy n = (x\ + \fdyi) n 
for n = 1, 2, 3 . . . 

Proof. First see that (x n , y n ) is a solution. We know that (x\,y\) is a solution 

x\ - dy 2 = I. 

Oi + \fdyi)(xi - Vdyi) = 1 

x n + \fdy„ = (x 1 + Vdyi) n 

taking conjugates, x n — Vdy n = {xi — Vdyi) n 

(x n + Vdy n )(x n - Vdy n ) = {x 1 + ^fdy 1 ) n {x 1 - Vdyi) n 

Xn - Vdy n = {{x\ + \fdyi)(xi - Vdyi)) n 
= l n = 1 



So (x n ,y n ) is indeed a solution. Why are these all the positive solutions? Sup- 
pose that (s, t) is a positive solution not of the form (x n , y n ) for any n. Then 
s + t\fd is a positive (> 1) real number. Not that {x n + Vdy„} is a sequence of 
positive real numbers which increase to infinity since 

I V 

x n + Vdy n = I xi +Vdyi I 

So pick n such that x n + y n Vd < s + t\fd < x n+ i + y n+ i\/d. Multiply the 
sequence of inequalities by x n — y n Vd (it's a positive real number because it 
equals x n - Vdy n = + ]/j y )• We see 

1 = (x n - Vdy n )(x n + Vdy n ) 

< (x n - Vdy n )(s + tVd) 

s v ' 

a+b^/d. a,&EZ 

< (x n - Vdy n )(x n+1 + Vdy n+1 ) 
^(x 1 -Vd yi ) n (x 1 + Vdy 1 ) n+1 
= (x 1 + Vdyi) 
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and so 

1 < a + bVd < xi + Vdyi 
We'll see a, b e N, then it will contradict minimality of {x\ , yi) 

a — bVd = = > 1 and a + bvQ > 1, so < a — bVd< 1 

a + bVd 

Adding 1 + < 2a gives a > \ > which means that a > 1. Also, 6 > ^ > 
=>■ (a, 6) is a positive integer solution. Why is a + 6\/d a solution? 

(a + bVd)(a - bVd) = (x n - Vdy n )(s + tVd)(x n + Vdy n )(s - tVd) 
a 2 -b 2 d=(xl-dy 2 n )( S 2 -dt 2 ) = l 



P-B equation is quite useful in many diophantine equations. 

Eg. Putnam asked, can we find infinitely many triples of consecutive integers, 
each of which is a sum of 2 squares? 

Yes. Suppose we choose n— 1, n, n+l where we set n = x 2 =>• n = x 2 +0 2 , n+l = 
x 2 + l 2 , x 2 - 1 = n - 1 = sum of 2 squares y 2 + y 2 , so we need to find infinitely 
many (x, y) such that x 2 — 2y 2 = 1. P-B, so ok. 

Proposition 79. Let N e Z, d e N, d ^ □. /fx 2 - dy 2 = N has one solution, it has 
infinitely many. 

Proof. Let (xi,y\) be a solution, so (xi + Vdyi){x\ — Vdyi) = N. Let (s n ,t n ) be 
infinitely many solutions to x 2 — dy 2 = 1 =>■ (s„ + Vdt n )(s n — Vdt n ) = 1. Then 
if we let x„ + v^y™ = + v^yijlsn + Vdt n ) it's easy to see x 2 — dy 2 = and 
that these are all distinct. So we get infinitely many solutions. ■ 

Eg. Prove that n 2 + (n + l) 2 is a perfect square for infinitely many values of n. 
Proof. 

n 2 + n 2 + 2n + 1 = 2n 2 + 2n + 1 = m 2 
An 2 + An + 2 = 2m 2 
(2n + l) 2 + 1 = 2m 2 

Let I be 2n + 1 =>■ get a solution of Z 2 + 1 = 2m 2 . Conversely, if I 2 + 1 = 2m 2 then 
I is odd, son — is an integer, and m 2 = n 2 + (n + l) 2 . (Just want to show 
that I 2 - 2m 2 = -1 has infinitely many solutions. We know it has an obvious 
solution (I, m) = (1,1) =>■ it has infinitely many.) ■ 
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Theorem 80 (Four Squares Theorem). Every non-negative integer is a sum of 4 
integer squares. 

Proof. Just like how we use complex numbers in the proof of the two squares 
theorem to establish that (a 2 + b 2 )(c 2 + d 2 ) = (ac - bd) 2 + (ad + be) 2 , we'll use 
quaternions now 

Q = {a + bi + cj + dk : a, b, c, d, € K} 
i, j, k are "imaginary" where 

i 2 = j 2 = k 2 = ijk = — 1 



ij = k ji = —k 

jk = i kj = —i 

ki = j ik = —j 

Multiplication in Q is non-commutative (but associative - x 1 (x 2 x 3 ) — (x\X2)Xz, 
etc). Addition is component- wise. If z = a + bj + cj + dk, define conjugate 

z = a — bi — cj — dk. Norm is ||z| = zz = a 2 + b 2 + c 2 + d 2 . 

Note that z~w-w-z.lt suffices to check things like 

-k = k = i] = j ■% = 

so 

\\zw\\ = zwzw = zwwz — z(ww)z — (zz)(ww) = ||z||||w| 

So (a 2 + b 2 + c 2 + d 2 )(e 2 + f 2 +g 2 + h 2 ) = (ae - bf - eg - dh) 2 + 3 other similar 
terms =>■ product of sum of 4 squares is a sum of 4 squares. So enough to show 
n is a sum of 4 squares for the case that n = 0, 1, or prime. = 2 + 2 + 2 + 2 , 
1 = l 2 + 2 + 2 + 2 , 2 = l 2 + l 2 + 2 + 2 , so enough to show that any odd 
prime p is a sum of 4 squares. 

Lemma 81. There's a positive integer m < p such that mp is a sum of 4 squares. 

Proof. Recall that if p is an odd prime then x 2 +y 2 + 1 = mod p has a solution 
(by pigeonhole principle). Let's suppose that we've produced x, y mod p so 

M,M<f.So 

x 2 + y 2 + 1< 2 (|) 2 + 1 = V - + 1< p 2 
So x 2 + y 2 + l 2 + 2 = mp for some < m < p. □ 

Let to be the smallest positive integer such that mp is a sum of 4 squares. We've 
showed m < p. If to = 1 done. So assume to > 1 and we'll get a contradiction 
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by producing a smaller value of m. If m is even, mp = x 2 + y 2 + z 2 + w 2 is 
even, so the number of odd elements of {x, y, z, w} is even. We can pair these 
up, say, as {x, y} and {z, w} such that x and y have same parity and z, w have 
same parity, so 

x + y x — y z + w z — w „ 

— -, — , e z 

2 ' 2 ' 2 ' 2 
x + j/^ 2 ^ / £ — y \ 2 ^ / z + ra\ 2 ^ ^ z — w ^ 2 .t 2 + y 2 + z 2 + w 2 



= ( ~2~) ^ (decreasing m) 



So suppose Z is some prime dividing m, necessarily odd. Write m = gl, so x 2 + 

y 2 + z 2 + w 2 = glp = mod I. Note I < p because l\m, m < p. Reduce x, y, z, w 
to x', y 1 , z', w' mod I (ie., x = x' mod I, etc.) such that \x'\, \y'\, \z'\, \w'\ < \. If 
x', y', z', w' are all 0, then x, y, z, w are also mod I, and 

(ir+(?r+(i) 2 +(M=f 

and j < gl = m reducing m. So we may assume x',y',z', w' are not all 0. 

x + yi + zj + wk = x' + y'i + z'j + w'k mod I 



Let 



Then 



p = x + yi + zj + wk 
a = x 1 + y'i + z'i + w'k 



WW = IMI 

= x' 2 + y' 2 + z' 2 + w' 2 

" v ' 

positive 

= x 2 + y 2 + z 2 + w 2 mod / 
= mod I 



So it's a multiple of I, say hi. Since \y\, \z\, \w\ < |, 



x 2 +y 2 +z 2 +w 2 < 4 ^J =l 2 
So < h< I. 

Also pa = pp mod I = x 2 + y 2 + z 2 + w 2 = mod I, so the components of 



2' 

I' " 
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quaternion pa are all divisible by I. Let /3 = e f. 

pa 

1 



= Ux 2 + y 2 + z 2 + w 2 )(x' 2 + v' 2 + z' 2 + w' 2 ) 

l Z v ^ ' v v ' 



hi 



(gip)(hi) 
l 2 



Note that m = gl and <?/i < gl since ft. < Z, so we have a sum of 4 squares which 
is a smaller multiple of p. ■ 
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